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Abstract 

Let Y be an Ornstein-Uhlenbeck diffusion governed by an ergodic finite state 
Markov process X: dY t — —X(Xt)Y t dt + a(Xt)dB t , Yq given. Under ergodicity condi- 
tion, we get quantitative estimates for the long time behavior of Y. We also establish a 
trichotomy for the tail of the stationary distribution of Y: it can be heavy (only some 
moments are finite), exponential- like (only some exponential moments are finite) or 
Gaussian- like (its Laplace transform is bounded below and above by Gaussian ones). 
The critical moments are characterized by the parameters of the model. 
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1 Introduction and main results 

The aim of this paper is to draw a complete picture of the ergodicity of Ornstein-Uhlenbeck 
diffusions with Markov switching (characterization of the tails of the invariant measure 
and quantitative convergence to equilibrium). In particular we make more precise the 
results of [7JI3]. The so-called diffusion with Markov switching Y = (Yt) t>0 is defined as 
follows. 

The switching process X = (Xt) t>0 is a Markov process on the finite state space 
E = {1, . . . , d} (with d ^ 2), of infinitesimal generator A = (A(x, x)) x ~ &E . Let us denote 
by a(x) the jump rate at state x S E and P = (P(x, x)) x xeE the transition matrix of the 
embedded chain. One has, for x ^ x in E, 

a(x) = —A(x,x) and P(x,x) — 



A(x, x) 



We assume that P is irreducible recurrent. The process X is ergodic with a unique 
invariant probability measure denoted by \i. See PHI for details. Let = a(X u , < u < 
t). Moreover, let E x denote the expectation with respect to the law ¥ x of X knowing that 
X = x. 

Let B = (Bt) t> Q be a standard Brownian motion on R and Yq a real-valued random 
variable such that B, Yq and X are independent. Conditionnally to X, the process Y = 
(Yt) t>0 is the real- valued diffusion process defined by: 

Y t = Y - f\{X u )Y u du + I a(X u ) dB u , (1) 
Jo Jo 
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where A and a are two functions from E to R and (0, oo) respectively. Of course, if A and 
a are constant, Y is just an Ornstein-Uhlenbeck process with attractive (A > 0), neutral 
(A = 0) or repulsive coefficient (A < 0). One has to notice that Equation (P) has an 
"explicit" solution: 



Y t = Y exp (- J \{X U ) dv)j + J exp (- J \{X V ) dv^j a{X u ) dB u . 



(2) 



Remark 1.1. In others words, the full process (X,Y) is the Markov process on E x R 
associated to the infinitesimal generator A defined by: 

Af(x, y) = J2 A & v) - /(*, v)) + ^j-9i 2 f(x, y ) - x(x)d 2 f(x, y ). 

Previous works investigated the ergodicity of Y and some integrability properties for 
the invariant measure. For example, in [2], the multidimensional case is adressed together 
with the case of diffusion coefficients depending on Y. Stability results and sufficient 
conditions for the existence of moments are established under Lyapunov-type conditions. 

In [7] , it is proved that the Markov switching diffusion Y is ergodic if and only if 

A(*)Mx) > 0, (3) 

xeE 

that is if the process is attractive "in average" . Let us denote by v its invariant probability 
measure of Y . It is also shown in [7] that v admits a moment of order p if, for any x £ E, 
p\{x) + a(x) is positive and the spectral radius of the matrix 

Msj^'H**, (4) 

is smaller than 1. In the sequel p(M) stands for the spectral radius of a matrix M. 

In [3], the result is more precise: a dichotomy is exhibited between heavy and light 
tails for v. Let us define 

A = minA(x) and A = maxA(x). (5) 

Theorem 1.2 (de Saporta-Yao [4J). Under Assumption (|3|) ; the following dichotomy 
holds: 

1. if A < 0, then there exists C > such that 

t K u((t, +oo)) ► C, 

t— >+oo 

where k is the unique p G (0, min {— a(x)/X(x), A(x) < 0}) such that the spectral 
radius of M p is equal to 1; 

2. if A ^ 0, then v has moments of all order. 
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Remark 1.3. Note that the constant k does not depend on the parameters {o~(x)) x£E , and 
that Point 1. from previous theorem implies that, for A < 0, the p th moment of v is finite 
if and only if p < k. 

The main idea of the proofs in [1] and /^/ is to study the discrete time Markov chain 
(Xg n , Y$ n ) n> Q for any 5 > with renewal theory and then to let 5 goes to 0. 

The main goal of the present paper is to show that there are three (and not only two) 
different behaviors for the tails of v. 

Let us gather below several useful notations. 

Notations 1.4. Let us define for the diffusion coefficients 

a 2 = miner 2 (x) and a 2 = max a 2 (x). (6) 
~ xeE xeE 

We denote by A p the matrix A — pA where A is the diagonal matrix with diagonal 
(A(l), . . . , A(ci)) and associate to A v the quantity 

rj p := — max Re 7. (7) 

7eSpcc(A p ) 

When A ^ 0, the set E is the union of 

M = {x e E, A(x) > 0} and N = {x G E, X(x) = 0}. (8) 
Let us then define 

o-(x) 2 - 

Mx) = -±J- and f3 = max/3(x), (9) 
2a(x) xeN 

and, for any v such that v 2 < (3 1 , the matrix 

We are now able to state our main result. 

Theorem 1.5. Let us define 

k = sup{p ^0, rj p > 0} G (0, +00]. 

Then r) p is continuous, positive on the set (0, k) and negative on (k, +00). Under Assump- 
tion ([3|) ; the following trichotomy holds: 

1. if A < then < k ^ min{— a(x)/X(x), X(x) < 0} ; and the p th moment of u is 
finite if and only if p < k. 

ifX = 0, then k is infinite and the domain of the Laplace transform of u is (—v c ,v c ) 
where 

v c = sup{v>0, p(Pi JV ))<l}; (11) 
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3. if A > 0, then k is infinite and v has a Gaussian-like Laplace transform: for any 
v € R, 



exp 



a 2 v 2 



4A 



=r- < e vy v{dy) < exp 



4A 



Moreover, its tail looks like the one of the Gaussian law with variance a/2 where 



a = max ie£ 



o~(x) 2 1 'A(x) since y i— > 



is v-integrable if and only if 5 < 1/a. 

and 3. as the 



Remark 1.6. In the sequel we will respectively refer to Points 1. 2. 
polynomial, exponential-like and Gaussian-like cases. 

The first point of this theorem is a reformulation of the first point of Theorem 11.21 by 
de Saporta and Yao. We can in particular check that our characterization of k in Theorem 
11.51 is equivalent to the one given by de Saporta and Yao in Point 1. of Theorem 11.21 (see 
Remark I4.3P . We provide a direct and simple proof of this result based on Ito formula 
and some basic results on finite Markov chains. The proof of Points 2. relies on precise 
estimates on the Laplace transform of Yt that can be derived from a discrete time model 
already studied in [6j [JJ . 

It is straightforward from ([2]) that, for any measure ttq on E x R, the Laplace transform 
L t of Y t is 

,2 



LAv) 



exp f«Y e- + J a(X, 



sj 2 e 



-2f*X(X r )dr ds 



(12) 



The estimate of the Laplace transform in the Gaussian-like case (Point 3.) 
easily deduced from this explicit expression. Assuming that Yq = 0, we get from 



is hence 
Uh that 



LAv) < E 



exp 



1. a 2 e~ 2 ^ dr ds 



< exp 1 



-2At 



4A 



which gives the upper bound as t goes to infinity. The lower bound follows from a sym- 
metric argument. 

The proofs of Point 2. and of the second part of Point 3. are more delicate (and 
interesting). For the exponential case, we first get the critical exponential moment for the 
process Y observed at the hitting times of the subset M defined in ([8]). Then we show 
that the full process has the same critical exponent. 

At the end of the paper we focus on the convergence of the law of Yt to the invariant 
measure v. We get an explicit exponential bound for the Wasserstein distance of order p 
for any p < k. Classically, let p ^ 1 and V v be the set of the probability measures on R 
with a finite p th moment. Define the Wasserstein distance W p on V v as follows: for any p 
and p in V p , 

i/p 



W p (p, p) = ( inf | J \y - y\ p n(dy, dy) 

where the infimimum is taken among all the probability measures tt on R 2 with marginals 
p and p. It is well-known that (V p , W p ) is a complete metric space (see [TT]). 

The strategy is to couple two processes (X,Y) and (X,Y) in such a way that the 
Wasserstein distance between C(Y t ) and C(Y t ) goes to zero as t goes to infinity. This 
requires to couple the initial conditions and the dynamics (of both the Markov chains and 
the diffusion part). When Xq and Xq have the same law, the coupling is trivial: we choose 
X = X and the same driving Brownian motion. 
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Theorem 1.7. Let p < k. Assume that Xq and Xq have the same law. Let Y and Y 
be solutions of ([T]) associated to (X t ) and (X t ) and assume that Yq and Yq have finite 
moment of order p. Then there exists C(p) such that 

wJc(Y t ),C(Y t )Y < C(p)e~^ t wJc(Y ),C(Y )'' P 



where r/ p is given by ([7]) . 

If Xq and Xq do not have the same law, one first has to make the Markov chains X 
and X stick together and then to use Theorem 1 1.71 This provides a rather intricate bound 
which is given for convenience in Section [5j 

The paper is organised as follows. In Section[2]we complete the proof for the Gaussian- 
like case of Theorem 11.51 The exponential-like case is studied in Section [3j Since the 
critical exponential moment is not explicit in the general case, we give also the explicit 
computation of the Laplace transform of v when E is reduced to {1,2}. In Section d] we 
establish a uniform bound for the p th moment of (Yt) t for any p < k and the first point of 
Theorem 11.51 as a corollary. We finally provide the proof of Theorem 11.71 and its extension 
to general initial conditions in Section 

2 Gaussian moments for the switched diffusion 

This section is dedicated to the proof of the second part of Point 3. of Theorem 11.51 
Proof of Point 3. of Theorem \1.5l Let us denote by 

a(x) = , , . for x £ E and a = maxa(x) < +oo. 
\(x) xeE 

For any 5 £ (0, 1/a), Ito's formula ensures that 

de SY ? = (-2\{X t )5Y t 2 + (25 2 Y t 2 + 5)a{X t f)e 5Y " dt + dM t 

where (Mt) t is a martingale. For any x £ E and y G R, 

2(-A(x) + 5a(x) 2 )y 2 + a(x) 2 ^ -2A(x)(l - 5a)y 2 + aX(x) 

^ -2A(1 - 5a)y 2 + aX, 

since 5a < 1. Moreover, for any a > 0, there exists b > such that, for any y G R, 

-2X5(1 - a5)y 2 t + Xa5 ^ -a + be- 5y2 , 

thus 

dt 



±E(e SY t 2 )^-aE(e 5Y ?)+b. 



As a consequence, sup^gE^e^ j is finite as soon as Efe a *° j is finite and 5a < 1. 

On the other hand, assume (without loss of generality) that a(l) = a. Choose (Xq, Yq) 
with law V (the invariant measure of (X, Y)). For any t > 0, we have 



E e SY o)=E e 5Y n >E 



l{X =i} E i,yo( 1 {r 1 >t}e 
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where T\ is the first jump time of X. On the set {T\ > t}, 

Y t = y e" A(1) * + N t 

where Nt is a centered Gaussian random variable with variance a(l)(l — e~ W*)/2 which 
is independent of Yq and T\. Thus, reminding that T\ ~ £(a(l)), we get 

Ei,yo(l { T 1>4} ^ 2 ) = e-( 1 )'E(e^- A(1)t ^) 2 ). 
Since a i— > E^e <5 ( a+Art * >2 ^ is even and convex, it reaches its minimum at a = and 

' if <5a(l)(l - e - 2A «*) < 1, 



+oo otherwise. 



As a consequence, if 5 > l/a(l), E^e 5y * j is bounded below by a function of t which is 
infinite for t large enoug h. Thus, E(V y * 2 ) is infinite too. □ 

3 Exponential moments for the switched diffusion 

This section is dedicated to the proof of Point 2. in Theorem 11.51 We assume in the 
sequel that A = 0. If (Xt) t>0 is a two-states Markov process then one can use (|12|) to 
compute explicitely the Laplace transform of the invariant measure v. This is a warm-up 
for the general case, and gives a more explicit formula for the critical exponential moment, 
whereas it will come from an abstract spectral criterion in the general case. 

3.1 The explicit expression for the two-states case 

In this subsection we assume that E = {1,2} and that A = 0. Let us start with a 
straightforward computation which suggests that the Laplace transform of the invariant 
measure of Y is infinite outside a bounded interval. 

Remark 3.1. IfT is an exponential random variable with parameter a and B is a standard 
Brownian motion on M. (with T and B independent) then, 



E(e vaBr )= / E(e vaBt )ae- at dt = 

Jo Jo 



e M/2 ae -at dt 



2a — a 2 v 2 



In other words, the law of gBt is a ( symmetric ) Laplace law. When X spends an expo- 
nential time in x G E with \{x) = 0, Y behaves like a{x)B. 

Theorem 3.2 (The two-states degenerate case). Assume that E = {1,2}, A(l) = A > 
and A(2) = 0. Then, for any v such that v 2 < 1/(3(2) (see ([9]) for the definition of (5), 

m = r e»^) = ( i -^m (r-La) 1+ " ,1,/A ^ cm) . 



l-P(2)v 2 J\l-p(2)v 2 J 4 ^ 
Ifv 2 ^ 1/(3(2), L(v) is infinite. 
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Proof. Since E = {1,2}, X is symmetric with respect to fx which is given by /x(l) = 
a(2)/(a(l) + a(2)). Let us denote by Lt the Laplace transform of Yt when Yq = and X 
is stationnary i.e. C(Xq) = fi. From Equation ([12]) . one has for any wgR, 



L t (v) = E, 



E, 



exp 



,2 /•* 



PQV 2 ^"^ 



exp j\(X s ) 2 e- 2 foMx r )dr ds 



since is reversible. By monotone convergence, we get that, for any veR, 

,,2 



L(v) = E^ 



exp J™a(X s ) 2 e- 2 SoHXr)dr d \ 



e [i,+oo], 



where L is the Laplace transform of v. 

Let us introduce two auxilliary functions: for x = 1,2, 



L x {v) = E a 



exp f^- J\(X s ) 2 e- 2 foKXr)dr ds \ 



It is clear that 



L(v)=fi(l)L 1 (v)+ f i(2)L 2 (v). 
Moreover, if for any t ^ 0, J-% = o~(X s , ^ s ^ t) and T is the first jump time of X, then 

,2 /-co 



La.(u) = E x 
= E X 



E x < exp 



v 



2 

2 /-T 



a(X s ) 2 e- 2 /oW* ds 



■Ft 



ew(jf °(X s ) 2 e- 2 fo^r) d r ds \ Ex 



where 



E x t = EJ exp 



,,2 />oo 



<7(X s ) 2 e- 2 J"o A P^d S 



For any s E [0, T[, X s = x and then 



a(X s ) 2 e' 2 ^ x ^ dr ds = a(x) 21 ' 



2A(x) ' 

with the convention (1 — e _0xT )/0 = T. Similarly, for t ^ T, 

^(X s ) 2 e- 2 i'o^)* ds = e -2A( a; )T / CT (x s ) 2 e- 2 ^W*d s 



The Markov property implies 

,2 /"oo 



E, 



Thus, 



eXP (T j^< X s) 2 ^ 2S ^ Xr)dr ds) 
v 2 a(x) 2 {l-e- 2X ^ T ) 



L x (v)=E 



exp 



4A(i) 



Xn= X 
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More precisely, 



Li(v) =Ei 



cxp 



v 2 a(l)\l-e 
4A 



-2AT N 



(„-»■) 



and 



L 2 (^)=E 2 r e ^ 2 ) 2T / 2 L 1 (t;)' 



2a(2) 



2a(2) - a(2) 2 v 2 

+ CX5 



if cr(2)V < 2o(2), 
otherwise. 



Using /3(2) = o-(2) 2 /2a(2), one easily gets that L\ satisfies the following equation: for any 
v 2 < 1/(3(2), 



1 f°° 

Li{v) = T^m?h exp 



1 - f3(2)v 2 Jo 



exp 



a{l) 2 v 2 (l-e- 2Xt ) 
4A 

\2„,2/ , i „,2> 



Li(t;e~ A *)a(l)e- a(1) *dt 



4A / A 



With x = uv, 

a(l)/A 



Ll(u) = _J__ Qy e -(l) 2 - 2 /(4A)^ e ^(l)^V(4A)^) x a(l)/A-l Li(a;)dx> 

Deriving this relation provides 



T> (v) - ( mV a(1) 4- U 4- 1 a W ^ r f„\ 
LlW " U"/3(2> 2 ~ "Aw" 2A + l-/3(2)^^rJ Ll W- 

Then Li is solution of 

^(1)^ /3(2)(l + a(l)/A)^ r 
LlW = (^A- + 1-/?(2V ) L1{V) 

which leads to 

/ 1 \ l+a(l)/A 

L 1 W=e CT ( 1 ) 2 ^/(4A)( 1 A 

1() \l-p(2)v 2 ) 

since -Li(O) = 1. Since L 2 is a function of L\ we get 

L(u) = e *(D^/(4A) f i-Ki)my \ ( l x 1+a(1)/A 



l-(3{2)v 2 J\l-P(2)v 2 i 

□ 
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3.2 The exponential-like case 

In this subsection we provide the proof of Point 2. (A = 0) of Theorem [T3J We first recall 
that, in this case, we split the state space E of the switching process X in two subsets M 
and N defined in ([5]). We denote also by F the points of M that can be reached in one 
step from N: 



Assume for simplicity that Xq E M and define by induction the sequence of times (T n ) n -^ Q 
by To = and, for n ^ 0, 

T 2n+1 = M{t > T 2n , X t £N}, and T 2n+2 = M{t > T 2n+1 , X t eM}. 

When X is in M , Y looks like a Ornstein-Uhlenbeck process (with variable but attractive 
drift) while it looks like a Brownian motion (with variable but bounded below and above 
variance) when X is in N. Thus, heuristically the process Y might be larger after a sojourn 
of X in N than in M. 

Let us notice that for x E N, 



and X x is the process X starting at x and T is the first hitting time of M. Our strategy 
is to determine the domain of the Laplace transform of I x and then to establish that is 
also the one of the process Y at the entrance times of X into the set M i.e at the times 
(T 2n ) n>0 . We will then extend the result to the full process (X,Y). 

Proposition 3.3. Under previous assumptions, for any v 2 < f3 1 , the two following 
conditions are equivalent: 

1. for any x E N, E(e vIx ) < +oo; 

2. p(Pv N ^) < 1, where Pv N ^ is defined in Equation (|1U|) . 

Proof. Let xq,x±, . . . ,x n _i be in N. We denote by (Z n ) n the embedded chain of X. On 
the set H = {Z = xq, . . . , Z n _i = x n -i, Z n E M}, 





where the 



£(r(xj)) 
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One just computes 

E(e wJ -o ) = Y, E ( eVlx ° \Zi = x 1 ,...,Z n _ 1 = x n - u Z n eM): 



x 1 ,...,x n - 1 eN 



x P X0 (^i = xi, . . .,Z n -i = x n -i,Z n G M)) 
P(x ,^i) P(x n _ 2 ,x n _i) P(x n -i,M) 



. 0(x o )v 2 1 - /?(x„_ 2 )u 2 1 - ^(x„-i)« 2 

Xl,...,Xn_ieJV 

% n-l 



n>l 



for y(ar) = 1 _^^ t)i P(%, M). Notice that tp is well-defined since v 2 < 1/(3. Moreover it is 
positive because X is irreducible recurrent, so, for any xq G N there exists a path that 
leads to M. 



If p(P„ (JV) ) < 1, then 



hmsup |^(pW)"-V| Vn < hmsup ||(^) B || 1/B < 1 , 

n— >+oo n^+oo 

hence the series is convergent. 

If p„ := p(Pv N ^) ^ 1, by Perron- Frobenius theorem, there exists a probability measure 
vq with some positive coefficients such that vqP^ = p v vo, which implies that 

E V0 (e vL )=u (<P)Y,Pv~ 1 = +oo > 

since ip is positive. □ 

Remark 3.4. When X is irreducible in restriction to N (i.e. the matrices Py N ^ are 
irreducible for any v), then K(e vIx ) = +oo for all x G N as soon as p(Pv N ^) ^ 1. If 
this it not the case, the previous proposition just ensures that when p(Pv N ^) ^ 1, then 
E(e vIx ) = +oo for some x G N. Moreover, for any x,x' G N such that P(x,x') is positive 
then E(e vI *') = +oo implies E(e vIx ) = +oo. 

We now introduce the sub-process made of the positions of (X, Y) at the successive 
hitting times of M. 

Proposition 3.5. For any n ^ 0, let us define 

U n = X T2n and V n = Y T2n . 
The process (U, V) is a Markov chain on F x R. More precisely, 

V n+1 = M n (U n )V n + Q n (U n ), 
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where the sequence of random vectors ((M n (x),Q n (x)) xeF ) is i.i.d., and independent of 
(U n ), with law given by 

M n (x) = exp (- f 1 \(X'*)dr 

Qn{x) 





\(Xf)exp \- £ \{X*)drJ dB s + j \{X x s )dB s . 

For any v < v c where v c = sup (v, p(Pv N ^) < l\, we have 

supEfe^ 1 ) < +00. 
Moreover, if v ^ v c , this supremum is infinite. 

Proof. The fact that (U, V) is a recurrent Markov chain is a straightforward application 
of the Markov property for X. 

Let us introduce M n = max x€ p M n (x) and Q n = max x€ p \Q n (x)\. The random vari- 
ables ((M n , Q n )) n>0 are i.i.d. Define the sequence (V n ) n>0 by 

V = \V \ and V n+ i = M n V n + Q n for n ^ 1. 

The domain of the Laplace transforms of (V n ) n>0 is known thanks to the exhaustive study 
PQ. Since P(Q n = 0) < 1, P(0 < M n < 1) = 1 and for any c£l, P(Q n + M n c = c) < 1, 
[TJ Theorem 1.6] ensures in particular that (Eexp (yVn)) n is uniformly bounded as soon 
as the Laplace transform L-q of Q is finite. At last, for any v ^ 0, 



supE(> IQ{:c)l ) <E(> Q ) = E^supe 1,|Q(:c) ^ sC J^K 6 " 



Thus Lq(v) is finite if and only if E(e 1, ' < ^^^) is finite for any x 6 F. Since ^ V n for 
all n ^ 0, then 

supE(V |Kl1 ) < +00 

as soon as L-q{v) is finite. 

On the other hand, choose v such that there exists xq £ i 7 such that E(e 1, ' < ^^ ^) is 
infinite. Then, for any re ^ 0, 



l {f7n=!Eo}e -l^l)E( e v l^^)l). 



> E 



The recurrence of U ensures that {n ^ 0, E(e v \ v "\) = +00} is infinite. 
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The last point is to show that Lq(v) is finite if and only if v < v c where v c is defined by 
(jlip . For any x G F, the random variable Q n (x) is symmetric and its Laplace transform 
is finite as soon as, for any x G N, the Laplace transform of 



h= [ T *(X*)dB s 
Jo 



is finite, which is true for \v\ < v c . Indeed, we have for any v 

=Xt 1 



Proposition 13.31 ensures that, if |f | < v c then 

E (>Q»(^ <; C^E^exp jT V(X;f) 2 exp ^(X^dr ) r/.s- 

Denoting <tm = max xg M <r(x) and A M = min xg A/ X(x), one has 

—2 



E (>Qn(*)) ^C7(«)exp (^« 2 ) 



By the way, Lq is finite on (— oo,v c ). 

We assume now that v ^ u c . From Proposition 13.31 we know that, in this case, the set 
G = {x G N, E(e vI «) = +00} is non empty. Using the irreducibility of X and Remark 
13.41 one notices that there exists xq G F such that ¥(X^P G G) > 0. From this remark 
and (JHJ), one has E(e t, ^™( x °)) = +00 which conclude the proof. □ 

Let us now extend this result to the whole process Y. 
Theorem 3.6. For any v < v c where v c = sup|w, p(Pv N ^) < l\, we have 

supE(e 1,|y ' |N ) < +00. 

Moreover, if v ^ v c , then this supremum is infinite. 
Proof. Choose t > 0. We have 



E ( e ^l)=£E( e ^ll {T2 ^ t<T2n 



t.+2j 

n=0 



We write, for ^ v < v c , 

^(e vlY %T 2n <t<T 2n+2} ) = E(E(e^ yt ll {r2ri<t<T2 „ +2} |^ T2n V^f 
As in the proof of Proposition 13. 5\ 

^ vlYtl MT^t<T 2n+2 }\^T 2n V^f ) < C(v)e W (^ 2 )eH^nlE(l {T2na<T2 „ +2} |^ T2 „ V^f). 
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By the Markov property applied to X, 

=2 



(e^ll {T2n<t<w ) ^ C(v)exp (^^ E ( e ^d)p(r 2n ^ t < T2n+2) 



E 

Then, for < v < v , 



E(e v W) < C(^)exp ( ^-v 2 ] supE(>K n 



The generalisation of the case v ^ v c to the whole process is immediate. □ 

4 Polynomial moments for the switched diffusion 

We denote by A p the matrix A — pA where A is the diagonal matrix with diagonal 
(A(l), . . . , A(<f)) and associate to A p the quantity 

rjp := — max Re 7. 

jeSpec(Ap) 

The main goal of this section is to establish the equivalence between the positivity of rj p 
and the existence of a p th moment for the invariant measure v of Y . We will also give the 
proof of Point 1 of Theorem 11.51 

Using classical ideas from spectral theory, we first relate r] p with exponential functionals 
of A along the trajectories of X: 

Proposition 4.1. For any p > 0, there exist < C\{p) < C2(p) < +00 such that, for any 
initial probability measure ir on E, any t > 0, 

CiCpJe-** < Eir (exp (- j\\{X u ) du^j \ < C 2 (p) e -**. (15) 

Proof. Let us define, for any p > and t > 0, the matrix -A^t) by 



A (p,t)( x > = E x (exp (- J p\{X u ) du^j t{ Xt =x) 



On the one hand, one remarks that 

Ejr (exp (- J p\(X u ) di?J \ = ttA m 1 (16) 

where the coordinates of 1 are all equal to 1 and ir is a probability measure on E seen as 
a row vector. 

On the other hand, a simple application of the Feynman-Kac formula shows that 
A/p f) = e tAp . This fact relates the spectra of A p and A^ pt y In particular, p{A^ p ^) = e~ Vpt 
and, since all coefficients of At p ^\ are positive, we can apply the Perron- Frobenius Theorem 
to ensure that — r/ p is a simple eigenvalue of A p , all other eigenvalues having a strictly 
smaller real part. Let £ p < — r\ p be an upper bound for the real parts of these other 
eigenvalues. 
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We then define ir p (resp. ip p ) the left (resp. right) eigenvector associated to — rj p , with 
positive coefficients, normalized such that vr p (l) = 1 (resp. ir p ((p p ) = 1). Applying [5j 
Thm VII. 1.8], we get that for any t > 

e tA p = e -Vpt (ppirp + R p (t) i 

with ||i? p (t)||oo ^ Pp(t)e^ pt , Pp{t) being a polynomial of degree less than d. This gives 

7re*M = e-**(7r(^p) + e^vr J R p (i)l) 

hence 

e-**(7r(yjip) - P p (i)e*(*+^) < 7re M n < e -**(7r(^ p ) + P p (i)e*^ + % 

This estimate gives (jT5|) thanks to (JTSD and to the fact that P p (t)e t{:n P +^ tends to as i 
tends to infinity. □ 

Let us now study the function p i— ► r] p . 

Proposition 4.2. 

1. The function p i— > ry p is smooth and concave on M + . /is derivative at p = is eguaZ 
to 

^ \(x)/j,(x) > 0, 

and np/p tends to X as p goes to infinity. 

2. We have the following dichotomy: 

• */ A ^ 0, i/ien /or a// p > 0, ?7 P > 0, 

• A < 0, i/iere is k £ (0, min{— a(x)/X(x), X(x) < 0}) sitc/i i/iai r/ p > /or 
p < k and rjp < /or p > k. 

Proof. The smoothness of the functions r/ p , 7r p and <£> p are classical results of perturbation 
theory (see for example [D] chapter 2]). Since ir p A p = —r] p 7r p , tt p 1 = 1 and Al = 0, one 
has 

rjp = -K p A p \ = pir p Al = p ^ k p {x)\{x). (17) 

Differentiating this relation gives rj p = 7r p Al + p7riAl. In particular, t/q = /xAl = 

T,xeEK x )K x )^ since = M- 

We turn to the proof of the concavity of r/ p . We only have to remark that, for any 
t > and any x £ E, 

p » M t (a °(p) = j lo S E x (exp f-p J \{X U ) du 

is a convex function, as a log-Laplace transform (for example using Holder's inequality). 
r x \ 

But (|15|) implies that M t converges to — rj p , hence rj p is concave as a limit of concave 
functions. 
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Obviously, one has, for any t > and p > 0, M^ x \p) ^ — pA and r\ p is greater than pA. 
On the other hand, denoting by T the first jump time of (Xt), one has 

Ml x \p) ^logljexp f-p jf A(X u )d«W >t} ) 

> -pA(sc) + jlogP z (T > t) = -pX{x) - a(x). 
When t goes to infinity, one gets for any p > 

7/p ^ min(a(x) + pA(x)). (18) 

In particular, %,/p goes to A as p goes to infinity. 

The fact that, when A ^ 0, r/ p is always positive is clear from (|17p . 

When A < 0, for p small enough, r] p > since its derivative at p = is positive. But in 
this case, we can check that rj p < for p large enough. Equation (|18p implies that r/ p < 
as soon as p > mm xeE i A( a; )<o ~~ a(x)/X(x). This provides the upper bound for k. 

With the concavity of r] p , these considerations are sufficient to ensure that Tjp as a 
unique zero k, being positive before and negative after. □ 

Remark 4.3. The relation rj K = implies that (A — K,K)(p K = which can be rewritten 
as M K ip K = <p K (M K being the matrix defined in flU) This ensures that p(M K ) = 1 since 
M K is non-negative irreducible and (p K is positive. By the way our characterization of k in 
Theorem \1.5\ is equivalent to the one given by de Saporta and Yao in Point 1. of Theorem 

ma 

It is known from 0] that the invariant measure v of Y has p th finite moment if and 
only if p < k. Their proof is based on a time discretization of the process (X, Y) together 
with generic results on the ergodicity of discrete time Markov processes and renewal theory 
(see |3j). The previous propositions provide a direct and simple characterization of the 
critical moment of v. 

Proposition 4.4. For any p > such that rj p > (i.e. p < k), and any initial measure 
such that the second marginal has a p th finite moment, one has 

supK(\Y t \ p ) < +oo and / \y\ p v(dy) < +oo. 
t^o J 

On the other hand, for any p such that rj p ^ (i.e. p ^ k) and any initial condition, 

lim E(|Yj| p ) = +oo and / \y\ p u(dy) = +oo. 
t^oo J 

Proof. Let us assume that p ^ 2. If it is not the case, one has to replace the function 



V l— * \y\ P by the C 2 function y i— ► 2 . Choose T > 0. Ito's formula ensures that 
d\Y t \ p = (-pX(X t )\Y t \ p + ^l-^la(X t ) 2 \Y t r 2 ^ dt + pa(X t )Y t \Y t f- 2 dB t . 



(19) 
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Let us denote by a p the function defined on [0, T] by 



a 



(*)=E(|Y t n^) 



Taking the expectation of (fT9j) conditionnally to X leads to 

a' p (t) = - P X(X t )a p (t) + P{P ~ 1) - a 2 (X t )a p - 2 (t), 
since B and X are independent. For any e > 0, there exists c such that 

u p (t) ^ (-pA(X t ) + e)ap(t) + c. 

This implies that 

•/ 

One has to take the expectation and use (|15|) to get for any p > 2 such that t/ p > 

n\Yt\ P ) < C 2 {p)n\Yo\ P )e { -^ +e)t + c C 2 (p) f e -(-*+ £ ) u du. 

Jo 

If e < rjp then sup t>0 E(|lt| p ) is finite. 
If p = k, one has 

a' K (t) = - K \(X t )a K (t) + K{K ~ l) a\X t )a^ 2 {t). 

Then 

a K (t) = / e~ K ^MXu)du^ K _ i) a (x s fa^ 2 {s)ds + ^{\Y Q \ K )e- K ^ x{x ^ du 
Jo 

rt t 

^k(k-1)(j 2 / e- K ^ A(x " )du a K _ 2 (s)ds. 
JO 

As a consequence, using Proposition 14. 1 1 and the relation rj K = (see Proposition 14. 2p . 

ft 

E{\Y t \ K ) ^ k(k - l)a 2 J E(a K „ 2 {s)E(e- K ^ x{x - )du \F^ds 
^ k(k - lk 2 Ci(K) J E(\Y s \ K ~ 2 ^j da. 
From the first part of the proof, 

Urn e(|F s |*~ 2 ) = j \y\ K - 2 v(dy) > 0. 

By the way, 

lim E(|Y t | K ) = +oo, 

t— >oo 

and the K th moment of v is infinite. This is also true for the p th moment for any p > k. □ 
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5 Convergence to equilibrium for the switched diffusion 



Under the assumption that v has a finite p th moment, one can establish an exponential 
convergence of (X, Y) to its invariant measure in terms of mixed total variation (for X) 
and W p Wasserstein distance (for Y). 

Let us start with the easiest case, assuming that C(Xq) = C(Xq). 

Proof of Theorem \t-l\ Let y and y be two real numbers. We couple two trajectories of 
(X, Y) starting at (x, y) and (x, y) by choosing the same first components and the same 
Brownian motion to drive Y and Y . In other words, we compare (Xt, Yt) x ' v and {Xt, Yt) x,y 
where 

Xt = Xt, 

Y t = y- [ \(X U )Y U du+ [ a(X u ) dB u 
J o Jo 

Y t = y- [ \(X u )Y u du+ [ a(X u )dB u . 
Jo Jo 



Then, 
and 



d Y-Y 



-X(X t )(Y t -Y t )dt 



Yt-Yt 



\y - y\ 



pKXu) 



Y —Y 

1 u 1 U 



du. 



As a conclusion, (fT5|) ensures that 



E 



Yt-Yt 



E x I exp 



pX(X u ) du\\\y- y\ p ^ C 2 (p)e-^\y - y 



Then, for any coupling LT of C(Yq) and C(Yq), 

W p (£(Y t ),£(Yt)y ^ C 2 (p)e-**y \y - y\ p U(d(y , y)) . 

Taking the infimum over II provides the result. 
Let us turn to the general case. 



□ 



Theorem 5.1. Consider two processes (X, Y) and (X,Y) with respective initial laws tt 
and tt two probability measures on E x R such that the second marginal has a finite th 
moment with 8 < k (with k = +oo if A ^ 0). For any p < 9, we have 



W p (£{Y t ),C{Y t )) < C 2 (p)(l -pc^Moiey^exp 
where 



{l-p/6)j + r) p 



t) + Pc Wte-^\ 



Pc = ^2fM)(x) A/2 (x) 
x&E 



l-d TV [C(Xo),C(X ] 

vie 



M (e)^ e = 2p(su P E(|y/) +su P E(|y/)] 

W p = m a xW p (c(Y \X = x),C(Y \X = x)), 
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and 7 is such that 



d TY (C(X t ),C(X t )) ^ e-^d TY (jC(X ),C(X )). 



Remark 5.2. This estimate can be improved and simplified if A > 0. In this case, one 
can write instead of (|20p that 



E 



(x,y),(x,y) 



Y t -Y t 



thanks to the explicit expression ([2]) ofY. Since pX ^ r] p this leads to 



W p [c(Y t ),C(Y t )Y sC C(p)(l-p c )exp 



Proof of Theorem \5.1{ We have to consider the case Xq ^ Xq. Given x,x G E (with 
i/i) and y,y G R, we introduce the three independent processes (X f )f^o, (^t)t>o an< ^ 
(Bt) t> Q where the first one is a chain starting at x, the second one is a chain starting at 
x and the last one is a standard Brownian motion. The process X is defined as follows: 



X t 



X t if t ^ T, 
X< if i > T, 



where T = inf ji > 0, X< = Xt}- It is well known (since X is a finite irreducible continu- 
ous time Markov chain) that there exists 7 > such that 

sup P X(5 (T > t) < e~ 7t . 

Let us now define for any t ^ 0, 

Let us denote, for any p < k and y,y € 



!>^ dv a{X u )dB u , 
f>^ )dv a(X u )dB u . 



C(p,x,y) = supE x , y (\Y t \ p ) and C(p, x, y, x, y) = 2 p (C(p, x, y) + C(p, x, y)). 
Let a G (0, 1) and s be the conjugate of 9 /p. Theorem 11.71 ensures that 



E 



(x,y),(x,y) 



E (x,y),(x,y) 



Yt-Y t 



1 



{T^at} + l{T<at} 



< C(e,x,y,x,y) p/e e-^ at/s 



(20) 



+ E 



(x,y),(x,y) 



Y T -Y T P C 2 (p)e~ v ^ t - T h 
(p)C(6, x, y, x, yfl e ( e ^ at l s + e^ 1 " )*) . 



{T<at} 
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Optimizing over a in order to have ^ya/s = r) p (l — a) i.e. a - „ 



E 



leads to 

IVp 



(|y t - y t [) *S C 2 (p)C((9, x, y, x, yY' 6 exp ("^^*) 



(x,j/),(x,S) I 



Let us now turn to the case of general initial conditions. Let ttq and ttq be two 
probability measures on E x R such that the second marginal has a finite th moment. 
Let us start coupling the marginals (io and jlo on E. Define the coupling probability p c 

Pc = ^2^o(x) Afio(x), 

and D = {x € E, no(x) ^ p,o(x)}. We introduce the random variables U, V, W and Z 
such that for any x e E 



F(U = x) = 
p(y = x ) = 
P(W = x) = 



Mo 0*0 A /io(x) 





Pc 


// (ar) 


- £o(z) 


1 


-Pc 


jl (x) 


- Po(x) 


1 


-Pc 



l D (x), 
1dc(x) 



and P(Z = 1) = 1 — P(Z = 0) = p c , Z being independent of (U, V, W). We can now define 



U ifZ = l, 

y if z = o, 



Xa = 



U if Z = 1, 
if Z = 0. 



We check by a standard computation that the law of Xo (resp. Xq) is /io (resp. fio). 

Now, for any x £ E, let us introduce two random variables Yq and Yq 1 , independent 
of (U, V, W, Z) such that 



E 



^ =^(£(y |X = a;),/:(lo|Xo = x)) 9 . 



With this construction (Xq,Y^°) has law ttq and (Xo,!^ ) has law ttq. We consider the 
processes (X, Y) and (X, Y) with these initial conditions, the sticky Markov chains and 
the same Brownian motion. Thanks to the previous computations, we have 

E(|y - y\ p ) = E(|y - Yt \ p (i {Xo=Xo} + 1 w *o})) 
^E(l {Xo=jto} \Y x ° - Y*°\y-^ 

+ C 2 (p)E(t {x ^ Xo} C\9,X ,Y x ^X X* o r /e ) exp (- 1Vp 
On the one hand, we have 



■t . 



{X =X () } 



^PcW^, 



E 



|Xn — Xi 



■)) 
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where W p = max l6 £ Wp^£(Yo|^o = x ), £(Yo\Xq = x)j . On the other hand, 

As a conclusion we get the following bound: 

wJc(Y t ),C(Y t )Y ^ C 2 (p)(l -p c )^M o (0)^e X p ( ^-t) +p 1 J s Wfe-^\ 

where 

M Q {ey' e = 2p(k(c(8,x ,y )) + ^{c{e,x ,%))) p/e . 

□ 
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